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Abstract 

We review the state of the art of f{R) theories of gravity (in their various formulations), 
which have been proposed as an explanation of the cosmic acceleration alternative to dark 
energy. The successes of f{R) gravity are discussed, together with the challenges imposed 
by minimal criteria for their viability. 
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1 Introduction 

The acceleration of the universe was discovered ten years ago using type la supernovae [1] and 
no definitive or truly satisfactory explanation of this phenomenon has been given yet. This 
discovery has important implications not only for cosmology, but also for fundamental physics. 
According to WMAP and the other experiments mapping anisotropics of the cosmic microwave 
background, if general relativity is the correct description of our universe, then approximately 
76% of its energy content is not dark or luminous matter, but is instead a mysterious form of 
dark energy, exotic, invisible, and unclustered. Three main classes of models for this cosmic 
acceleration have been proposed: 

1. a cosmological constant A 

2. dark energy 

3. modified gravity. 

Naively, a cosmological constant propelling the cosmic acceleration and eventually coming to 
dominate the universe, causing it to enter a de Sitter phase without return, seems the most 
obvious explanation. However, A brings with it the notorious cosmological constant problem 
and the coincidence problem. To admit that A is non-zero but still has the tiny value required to 
cause the current acceleration amounts to an extreme fine-tuning. It is not surprising, therefore, 
that most cosmologists reject this explanation, postulating that A is exactly zero for reasons yet 
to be discovered (the cosmological constant problem is put in a drawer for the time being), and 
that a different explanation is to be found for the cosmic acceleration. 

The second class of models, (mostly) within the context of general relativity, postulates the 
existence of a dark energy fluid with equation of state P ~ —p (where p and P are the energy 
density and pressure of the fluid, respectively), which comes to dominate late in the matter era. 
Dark energy could even be phantom energy with equation of state such that P < —p. Many 
dark energy models have been studied, none of which is totally convincing or free of flne-tuning 
problems, or can be demonstrated to be the "correct" one. 

A third possibility consists in dispensing entirely with the mysterious dark energy and mod- 
ifying gravity at the largest scales0 Here we focus on modifled gravity and, speciflcally, on the 
so-called f{R) gravity theoriesH 

f{R) or "modifled" gravity consists of infrared modiflcations of general relativity that become 
important only at low curvatures, late in the matter era. The Einstein-Hilbert actiorH Seh = 
^ / d^x,/^R + is modifled to 

S = ^ [ d^x^f{R) + S^"^"''^'^ , 
Ik J 

where f{R) is a non-linear function of its argument OH]. 

In principle, the metric tensor contains several degrees of freedom: tensor, vector, and scalar, 
massless or massive. In general relativity only the familiar massless spin 2 graviton propagates. 
When the Einstein-Hilbert action is modified, other degrees of freedom appear. The change 

^Historically, this approach was the correct one in explaining the precession of Mercury's perihelion: not due 
to an unseen mass but to Einstein's modification of Newtonian gravity. 

^We do not consider here braneworld models or other modifications of gravity that are sometimes also referred 
to as "modified gravity". 

^Here R is the Ricci curvature of the metric tensor gab with metric determinant g, k = SttG, G is Newton's 
constant, and we follow the notations of Ref. [2]. 



(1.1) 



i? — 2A — > f{R) in the action brings to hfe, in addition to the massless graviton, a massive scalar 
mode which can drive the cosmic acceleration, and will be discussed below. This is analogous 
to the inflaton field driving the accelerated expansion of the early universe, although at a much 
lower energy scale. 

If terms quadratic in the Ricci and Riemann tensor, and possibly other curvature invariants, 
are included in the gravitational Lag rangian, f i^R^ Rq})R^^ ^ Rqj^q^R^^'^'^ ^ •••); massive gravitons 
and vector degrees of freedom appear. In the following we restrict ourselves to f{R) theories of 
gravity in their various formulations, and we focus on their use as substitutes for dark energyQ 
for a more comprehensive discussion we refer the reader to |7j and for short introductions to 

mm- 

We adopt a conservative point of view and regard f{R) gravity more as a toy model than the 
correct theory of gravity, i.e., we consider these theories as a proof of principle that modifying 
gravity is a viable alternative to dark energy. However, we do not feel that one can claim 
that any of the f{R) models proposed thus far is the "correct" one, or has exceptional support 
from the observational data. While it is true that many f{R) models pass all the available 
experimental tests and fit the cosmological data, the same is true for many dark energy models, 
and it is currently impossible to use observational data to discriminate between most of them, 
and between dark energy and modified gravity models. 

Modifying gravity is risky: unwanted consequences may be violations of the experimental 
limits on the parametrized-post-Newtonian (PPN) parameters at terrestrial and Solar System 
scales [To], instabilities, ghosts and, as in any newly proposed theory, the Cauchy problem could 
be ill-posed. These aspects are discussed in the next sections. 

f{R) gravity has a long history: its origins can be loosely traced to Weyl's 1919 theory 
in which a term quadratic in the Weyl tensor was added to the Einstein-Hilbert Lagrangian 
|11| . Later, f{R) gravity received the attention of many authors, including Eddington, Bach, 
Lanczos, Schrodinger, and Buchdahl. In the 1960's and 1970's, it was found that quadratic 
corrections to Seh were necessary to improve the renormalizability of general relativity |12j . 
and in 1980 quadratic corrections were found to fuel inflation without the need for scalar fields 
|13j . Non-linear corrections are also motivated by string theories |14j . We refer the reader to 
|15j for an historical review. 

The prototype of f{R) gravity [3111] is the model 

/(i?) = i? - /iV^ , (1.2) 

where /i is a mass scale of the order of the present value of the Hubble parameter /x ~ Hq ~ 
10~^^ eV. Although ruled out by its weak-field limit [16] and by a violent instability [17], this 
model gives the idea underlying modified gravity: the 1/R correction is negligible in comparison 
with R at the high curvatures of the early universe, and kicks in only as ii ^ 0, late in the 
history of the universe. 

Many forms of the function f{R) are found in the literature: here we discuss only general, 
model-independent, features of f{R) gravity. 

2 The three versions of f{R) gravity 

Modified gravity comes in three versions: 
1) metric (or second order) formalism; 



f{R) gravity, sometimes with an explicit coupling of matter to -R has also been used as an alternative to 
galactic dark matter [6]. 



2) Palatini (or first order) formalism; and 

3) metric-affine gravity. 



2.1 Metric f{R) gravity 

In the metric formalism [3l 0] , the action is 



^metric 



1 

2k 



Variation with respect to the (inverse) metric tensor g*^^ yields the field equation 

f{R) 



f'{R)Ra 



■ 9ab = VaVb/'(i?) - gabDf'iR) +KTab 



(2.1) 



(2.2) 



where a prime denotes differentiation with respect to R. The first two terms on the right hand 
side introduce fourth order derivatives of the metric, hence the name "fourth order gravity" 
sometimes given to these theories. 
The trace of eq. ()2.2p yields 



3nf{R)+Rf{R)-2fiR) = KT 



(2.3) 



where T = T^a is the trace of the matter stress-energy tensor. This second order differential 
equation for f'{R) differs deeply from the trace of the Einstein equation R = —kT which, 
instead, relates algebraically the Ricci scalar to T. We already see that f'{R) is indeed a 
dynamical variable, the scalar degree of freedom contained in the theory. 

Formally, one can rewrite the field equation ()2.2|) in the form of an effective Einstein equation 

as 



where 



pieff) 
' ab 



Gab 

f{R) - Rf'iR) 



-{Tab + Tj^l"'^ 



gab + VaVbfiR) - gabOf'iR) 



(2.4) 



(2.5) 



is an effective stress-energy tensor containing geometric terms. Of course, as usual when adopting 
this procedure, T^^^^ does not satisfy any energy condition and the effective energy density is, 
in general, not positive-definite. As is clear from these equations, in f{R) gravity one can define 
an effective gravitational coupling Ge// = G/ f'{R) in a way analogous to what is done in scalar- 
tensor theories. Hence, f'{R) must be positive in order for the graviton to carry positive kinetic 
energy. 

In the spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) metric adopted as the 
kinematic description of our universe, 

ds^ = -df + a^{t) [dx"^ + dy'^ + dz^) , (2.6) 

the field equations of metric f{R) cosmology assume the form 

Rf'iR) - f{R) 



2H + 3i?2 = 



P 



(matter) 



+ 



- 3HRf"{R) 



(2.7) 



f{R)-RnRy 



pimatter) ^ f p\ ^ 2HRf"{R) + Rf" {R) 



(2.5 



where an overdot denotes differentiation with respect to the comoving time t. The corresponding 
phase space is a 2-dimensional curved manifold embedded in a 3-dimensional space and with a 
rather comphcated structure [I8j . 



2.2 Palatini f{R) gravity 

In the Palatini approach, both the metric gab and the connection T^^ are independent variables, 
i.e., the connection is not the metric connection oi gab- While in general relativity the metric and 
Palatini variations produce the same (Einstein) equations, this is no longer true for non-linear 
LagrangiansH 

Shortly after metric f{R) theories were proposed as alternatives to dark energy, also the 
Palatini version was adanced for the same purpose, originally in its f{R) = R—fi'^/R incarnation 
|20j . The Palatini action is 

Spalatim = J d^'x ^ f (ij) + S^"""^'') [ffafc,^'^"^^] • (2-9) 

There are two Ricci tensors: Rab, which is constructed using the metric connection of the 
(unique) physical metric gab, and Rab which is the Ricci tensor of the non-metric connection F^^. 
The latter gives rise to the scalar R = g'^^Rab- The matter part of the action does not depend 
explicitly from the connection T, but only from the metric and the matter fields, collectively 
denoted with ^/;(™). 

Variation of the Palatini action (j2.9p yields the field equation 



r{R)Rab - ^ gab = KTab. (2.10) 

Note the absence of second covariant derivatives of /', in contrast with eq. (j2.2p . Variation with 
respect to the independent connection produces the field equation 

Vd (^nV') - Vd (^/'(i?)/^'^) = , (2.11) 



where Vc denotes the covariant derivative associated to this non-metric connection T. The trace 
of eqs. ([XTUD and OTTD yields 

f'{R)R-2f{R)=KT (2.12) 

and 

vJv^f{R)g''')=0. (2.13) 



Eq. (|2.13|) tells us that Vc is the covariant derivative of the metric 

gab = f{R)gab (2.14) 

conformally related to gab- Note that eq. (j2.12p is an algebraic (or trascendental, depending on 
the form of the function /) and not a differential equation for f'{R): hence, this quantity is 
non-dynamical, contrary to metric f{R) gravity. This lack of dynamics has important conse- 
quences explored in the following sections. It is possible to eliminate completely the non- metric 



^The requirement that the Palatini and metric variations give the same field equations selects Lovelock gravity 
[19) . of which general relativity is a special case. 



connection from the field equations, which are then rewritten as 



Gab j,/ T^ab 
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(2.15) 



2.3 Metric-afflne f{R) gravity 

In metric- affine f{R) gravity [21j, also the matter part of the action 

Saffine = ^ J d\ ^ f [r) + S^^'^''^'^ [^af,, T^, V'^")] , (2.16) 

depends explicitly on the connection F, which is possibly non-symmetric. This leads to a torsion 
associated with matter, and to a modern revival of torsion theories. These were originally 
introduced within a non-cosmological context, with the spin of elementary particles coupling to 
the torsion. Metric-affine f{R) gravity has not yet been explored in great detail, especially with 
respect to its cosmological consequences. For these reasons, in the following we focus on metric 
and Palatini f{R) gravity. 

3 Equivalence of metric and Palatini f{R) gravities with Brans- 
Dicke theories 



If f"{R) 7^ 0, metric modified gravity is equivalent to an = Brans-Dicke theor5|j [22], while 
Palatini modified gravity is equivalent to an w = —3/2 one. This equivalence has been proposed 
and rediscovered, for particular theories or in general, many times over the years |23J. 

3.1 Metric formalism 

Assuming that f"{R) ^ and beginning with the action (II. ip . one introduces the auxiliary 
scalar field cj) = R and considers the action 

S = ^^j d^x^g [V(0)i? - Vm + S^^^''-^'^ , (3.1) 

where 

^(<A) = /'(<A) , V{ck) = ^f'icp) - /(-/>) . (3.2) 

The action (j3.ip trivially reduces to (jl.ip for metric f{R) gravity (j) = R. Vice- versa, the 
variation of (|3.ip with respect to g""^ gives 

1 / V \ K 

Gab = — [ ^aSbi^ - gabOl/j - — gab) + - Tab , (3.3) 



while varying with respect to (j) yields 



^^-^ = (^-'A)m = (3.4) 



''The general form of the Brans-Dicke action is Sbd ^ ^ ! d'^x \j)R - ^ - + 



and (f) = R under the assumption /" 7^ 0. Hence, the scalar field <j) = R is dynamical and 
satisfies the trace equation 

3/"(0)n0 + 3/"'(0) V^0Ve<^ + (Pfi^) -2f{^)=KT. (3.5) 

This scalar is massive: as discussed later, the analysis of small perturbations of de Sitter space 
allows one to compute explicitly its mass squared 

^ = - ^0) , (3.6) 

where a zero subscript denotes quantities evaluated at the constant curvature of the de Sitter 
background. It turns out to be more convenient to consider the scalar ip = f'{(p), which satisfies 

SDV' + 2U(ip) -i;^ = KT (3.7) 
dip 

with U{^) = V{(l){ip)) — f{(f){tjj)). It is clear, therefore, that the theory contains a scalar degree 
of freedom, and the action 

S= — I (fx ^ [tPR - U{i;)] + S^'"""^'') , (3.8) 
2k J 

is recognized as an = Brans-Dicke theory. This theory, called "massive dilaton gravity" 
was originally introduced in the 1970's in order to generate a Yukawa term in the Newtonian 
limit [23]. The assumption /" 7^ can be seen as the requirement that the change of variable 
R i^iR) be invertible. 



3.2 Palatini formalism 



In the Palatini case, the discussion of the equivalence with a Brans-Dicke theory proceeds in a 
way analogous to that of the metric formalism. One begins with the action (|2.9|) and introduces 
(j) = R and ij) = /'(</>)• Then, apart from a boundary term that can be neglected for classical 
purposes, the action is rewritten, in terms of the metric gab and of its Ricci tensor Rat, as 



S 



Palatini 



1 

2^ 



_|_ g(matter) 



(3.C 



where we used the fact that, since gah = '4' 9 ah-, the Ricci curvatures of gah and gah are related by 

(3.10) 

The action (j3.9p is recognized as a Brans-Dicke theory with Brans-Dicke parameter uj = —3/2. 



R = R+^\/^^l;\Jc^-lu^ . 
lip 2 



4 Criteria for viability 

In order for f{R) gravity to be successful, it is not sufficient that it serves the purpose for which 
it was introduced in the cosmological context, but it must also pass the tests imposed by Solar 
System and terrestrial experiments on relativistic gravity, and it must satisfy certain minimal 
criteria for viability. Overall, these are: 



possess the correct cosmological dynamics; 



• not suffer from instabilities and ghosts; 

• have the correct Newtonian and post-Newtonian limit; 

• give rise to cosmological perturbations compatible with the data from the cosmic microwave 
background and large scale structure surveys; and 

• have a well-posed Cauchy problem. 

The failure to satisfy even a single one of these criteria is taken as a statement that the 
theory is doomed. These viability criteria are examined in the following. 

4.1 Correct cosmological dynamics 

In the opinion of most cosmologists, in order to be acceptable a cosmological model must exhibit 
early inflation (or an alternative way to solve the horizon, flatness, and monopole problem to- 
gether with a mechanism to generate density perturbations), followed by a radiation-dominated 
era and a matter-dominated era, and then by the present accelerated epoch that f{R) theories 
were resurrected to explain. The future era is usually found to be an eternal de Sitter attractor 
phase, or a Big Rip singularity truncating the history of the universe at a finite time. 

Smooth transitions between different eras are required. It has been pointed out that the 
exit from the radiation era, in particular, may have problems in many models [25], a warning 
that care must be exerted in building f{R) cosmologies. Ultimately, exit from the radiation, or 
any era can be achieved. Take, for example, what we could name "designer f{R) gravity": one 
can prescribe a desired expansion history of the universe by a choice of the scale factor a{t) and 
then integrate an ODE that determines the function f{R) that produces a{t) [26]. In general, 
this function is not unique and assumes rather contrived forms (not the usual R — /U^*-""''^-'//?", 
or simple forms like that). 

4.2 Instabilities 

The prototype model in the discussion of instabilities is again the choice f{R) = R - ^i^/R 
with Ho ^ 10~33 eV. Shortly after it was proposed, this model was found to suffer from a 
catastrophic ( "Dolgov-Kawasaki" ) instability [l^. The stability analysis was later generalized 
to any metric f{R) theory [27J and the extension to even more general gravitational theories 
has been pursued [28]. One proceeds by parametrizing the deviations from general relativity as 




where e is a small positive constant with the dimensions of a mass squared and the function ip 
is dimensionless. The trace equation for the Ricci scalar R takes the form 




f{R) = R + e^{R) 






Next, one expands around a de Sitter background and writes the metric locally as 



9ab — Vab + hab 



(4.3) 



while the scalar degree of freedom R is expanded as 



R = -kT + Ri , 



(4.4) 



with Ri a perturbation. To first order, the trace equation yields the dynamical equation for Ri 



R^-VR^-^TR, + —^VT.VR, + —i^--^jR, = KT-KVT — . 

(4.5) 

The last term on the left hand side is dominated by the term in e ^ and gives the effective mass 
squared of Ri 

- ^ , (4.6) 

from which one deduces that the theory is 

• stable if f"{R) > 

• unstable if f"{R) < . 

The case of general relativity is excluded by the assumption f'^O, but the well-known stability 
in this case allows one to extend the stability criterion for metric f{R) gravity to be /" > 0. 

As an example, the prototype model f{R) = R — fJ-'^/R, which has /" < is unstable. The 
time scale for the onset of this instability is dictated by the smallness of the scale and is seen 
to correspond to ^ 10~^^ s [T7], making this an explosive instability. 

One can give a physical interpretation of this result as follows [29]: remembering that the 
effective gravitational coupling is Gg// = G/f'{R), if dGeff/dR = -f"G/{f')'^ > (which 
corresponds to /" < 0), then Ge// increases with R and a large curvature causes gravity to 
become stronger, which in turn causes a larger R, in a positive feedback mechanism driving the 
system away. If instead dGeff/dR < 0, then a negative feedback damps the increase in the 
gravitational coupling strength. 

Palatini f{R) gravity, by contrast, is described by second order field equations, the trace 
equation f'{R)R — 2f{R) = kT is not a differential equation but rather a non-dynamical alge- 
braic one and, therefore, there is no Dolgov-Kawasaki instability [30j. 

The previous analysis for metric f{R) gravity obtained with the local expansion (14. 3|) is 
necessarily limited to short wavelengths (compared to the curvature radius) , but can be extended 
to the longest wavelenghts [31j . This is necessarily more complicated because these modes suffer 
from the notorious gauge-dependence problems of cosmological perturbations and a covariant 
and gauge-invariant formalism is needed. We assume that the background space is de Sitter and 
consider the general action 



/ 



S = d X \/—g 



,4 ^ /(</., i?) ^(0) 



(4.7) 



which contains f{R) and scalar-tensor gravity, and mixtures of them. On a FLRW background, 
the field equations become 

H' = ^[^^' + ^^ + V-^Hf^ , (4.8) 
H=^(u:^' + f'-Hf') , (4.9) 



2f' 



(4.10) 



de Sitter space is a solution subject to the conditions 



GH^fo -fo + 2Vo = 0, f^ = 2^0' . (4.11) 

An analysis [31] using the covariant and gauge-invariant Bardeen-Ellis-Bruni-Hwang formalism 
[32] in the version given by Hwang [33] for alternative gravitational theories yields the stability 
condition of de Sitter space in metric f{R) gravity with respect to inhomogeneous perturbations 

^^°^'7„y°^° > , (4.12) 

JOJO 

which is obtained in the zero momentum limit. This condition coincides with the stability 
condition with respect to homogeneous perturbations [29]. 

At this point it is worth checking that the equivalence between metric f{R) gravity and an 
a; = Brans-Dicke theory holds also at the level of perturbations; previous doubts to this regard 
[34] [35] have now been dissipated. 

For the a; = Brans-Dicke theory, the stability condition of de Sitter space with respect to 
inhomogeneous perturbations is given again by eq. (j4.12p . while that for stability with respect 
to homogeneous perturbations is 

{fo? - 2hfl! > _ (413) 
/o 

This is again equivalent to (|4.12|) provided that stability against local perturbations, expressed 
by /o > 0, is assumed. Therefore, there is complete equivalence between metric f{R) gravity 
and iv = Brans-Dicke theory also at the level of cosmological perturbations. 

Going beyond the linear approximation, metric f{R) theories have been found to be suscep- 
tible to another, non-linear, instability, which makes it hard to build models of relativistic stars 
in strong gravity. In these situations, a singularity develops for large R, which was discovered in 
[36j . Although this problem needs further study, it seems that, in order to avoid this singularity 
requires some degree of fine-tuning. At present, this is probably the biggest challenge for metric 
f{R) theories. 



4.3 Ghosts 

Ghosts are massive states of negative norm which cause lack of unitarity and are common when 
trying to generalize Einstein's gravity. The good news here are that f{R) gravity is ghost-free. 
More general theories of the form / (^R, RabR"'^ , RabcdR"'^'^'^ ^ ■■■) ^ in general, contain ghost fields. 
A possible exception (under certain conditions [37]) is the case in which the extra terms appear 
in the Gauss-Bonnet combination Q = R^ - ARabR"-^ + RahcdR"'^"'^ , as in / = / (i?, g) In this 
case, the field equations are of second order and there are no ghosts [38 ] [39 l [30 ] . 



4.4 Weak-field limit (metric formalism) 

Early work on the weak-field limit of both metric and Palatini f{R) gravity was subject to 
errors and incompleteness (see [7] for detail^; a satisfactory treatment for the prototype model 
f{R) = R — n'^/R in the metric formalism was given in [42| and then generalized to arbitrary 
forms of the function f{R) in [43\ 116] . 



^The limit of f{R) gravity to general relativity has the character of a singular limit, and even the limit to 
general relativity of Brans-Dicke theory is not free from ambiguities [41) . 



In order to assess whether the hmits set on the PPN parameter 7 by the available Solar 
System experiments, one needs to find the weak-field solution of the field equations and com- 
pute this parameter. One considers a static, spherically symmetric, non-compact body which 
constitutes a perturbation of a background de Sitter universe. The line element is written as 

- [1 + 2^'(r) - Hlr"^] dt^ + [l + 2$(r) + i^^r^] dr"^ + r^dO^ (4.14) 

in Schwarzschild coordinates, where is the line element on the unit 2-sphere and ^' and 
$ are post-Newtonian potentials. These are of small amplitude, , << 1, and one 

considers small (non-cosmological) scales so that Nqv « 1, while expanding the Ricci scalar 
around the constant curvature of the background de Sitter space, R{r) = Rq + Ri. The PPN 
parameter 7 is given by 7 = — <I>(r)/^'(r) [lOj. Three assumptions are made |16j : 

1. f{R) is analytical at Rq; 

2. mr « 1, where m is the effective mass of the scalar degree of freedom of the theory, i.e., 
it is assumed that this scalar field is light and has a range larger than the size of the Solar 
System (we remind the reader that there are no experimental constraints of scalars with 
range m^^ < 0.2 mm). 

3. For the matter composing the spherical body, the pressure is negligible, P ~ 0, so that 
T = To + Ti ~ -p. 

The first and the last assumption are not stringent, but the second one is, as will be clear 
below. The trace equation p.3p yields the equation for the Ricci scalar perturbation 



V^Ri - m^Ri = , (4.15) 
-J Jo 



where 

= (4.16) 

is the effective mass squared of the scalar. Eq. ()4.16p coincides with the expression obtained in 
the gauge-invariant stability analysis of de Sitter space and in propagator calculations. 
If mr < < 1 , the solution of the linearized field equations is 

, , -kM 1 

= IJ-TT - > (4-17) 



6vr/o r 
kM 1 



/ , K.M 1 
Hr) = TTT^ - ■ (4.18) 



The PPN parameter sought for is, therefore. 



-^>(r) 1 

in gross violation of the (recently improved) experimental limit |44j 

I7 - 1| < 2.3 • 10"^ . (4.20) 



This result would be the end of metric f{R) gravity if the assumptions made in the calculation 
were satisfied. However, this is not the case for assumption 2): mr is not always less than unity 



due to the chameleon effect. This consists in the effective mass m depending on the curvature 
or, alternativeiy, the matter density of the environment. The scalar degree of freedom can be 
short-ranged (say m > 10~^ eV, corresponding to a range A < 0.2 mm) at Solar System densities 
and evade the experimental constraints, while being long-ranged at cosmological densities and 
thus being able to affect the cosmological dynamics [301 US]. Although at a first glance the 
chameleon effect could be seen as a contrived and fine-tuned mechanism, f{R) gravity is rather 
complicated and the effective range does indeed depend on the environment. The chameleon 
mechanism is well-known and accepted in quintessence models, in which it was discovered for 
the scalar field potential V{(j)) ~ !/(/> [l6]. Many forms of the function f{R) are known to exhibit 
the chameleon mechanism and pass the observational tests. For example, the model 



is compatible with the PPN limits if ~ lO"^'^ eV~ 10^^^ Hq [45] . It is obvious that a correction 
term ~ i?" with n < 1 to the Einstein-Hilbert Lagrangian R will come to dominate as i? ^ 0"*" 
(for example, for n = 1/2, > as i? — > 0). The model ()4.2ip is compatible with the 
experimental data but it could be essentially indistinguishable from a dark energy model. Hope 
of discriminating between dark energy and f{R) models, or between different modified gravities 
relies on the study of the growth of cosmological perturbations. 

4.5 Correct dynamics of cosmological perturbations 

The expansion history of the universe alone is not sufficient to discriminate between various 
models, but the growth of structures depends on the theory of gravity and has the potential to 
achieve this goal. Song, Hu, and Sawicki [47] assumed an expansion history a{t) typical of the 
ACDM model and found that vector and tensor modes are not affected by f{R) corrections to 
Einstein gravity, to lowest order, while scalar modes are. They also found the condition f"{R) > 
for the stability of scalar perturbations, in agreement with the arguments discussed above. 
The most interesting results are that f{R) corrections lower the large angle anisotropics of the 
cosmic microwave background and produce correlations between cosmic microwave background 
and galaxy surveys different from those of dark energy models. 

Overall, the study of structure formation in modified gravity is still work in progress, and 
often is performed within the context of specific models, some of which are already in trouble 
because they do not pass the weak- field limit or the stability constraints. A similar situation 
holds for all Palatini f{R) models, and for this reason, their weak-field limit and cosmological 
perturbations are not discussed here. 

4.6 The Cauchy problem 

A physical theory must have predictive power and, therefore, a well-posed initial value problem. 
General relativity satisfies this requirement for "reasonable" forms of matter [2]. The well- 
posedness of the Cauchy problem for vacuum f{R) gravity was briefly discussed for special 
metric models in earlier papers [38]. Thanks to the equivalence between f{R) gravity and 
scalar-tensor theory when f"{R) ^ 0, the Cauchy problem can be reduced to the analogous 
one for Brans-Dicke gravity with uj = 0, —3/2. That the initial value problem is well-posed was 
demonstrated for particular scalar-tensor theories in [49 1 BH] and a general analysis has only 
recently been performed [50[ I51j . This work, however, does not cover the a; = 0, —3/2 cases. 




(4.21) 



A system of 3 + 1 equations of motion is well-formulated if it can be written as a system of 
equations of only first order in both time and space derivatives. If the latter is cast in the full 
first order form 

dtu + M'ViU = S{u), (4.22) 

where u collectively denotes the fundamental variables hij,Kij, etc. introduced below, M* is 
called the characteristic matrix of the system, and S (u) describes source terms and contains only 
the fundamental variables but not their derivatives. The initial value formulation is well-posed if 
the system of PDEs is symmetric hyperbolic (i.e., the matrices M* are symmetric) and strongly 
hyperbolic if SjM' has a real set of eigenvalues and a complete set of eigenvectors for any 1-form 
Si, and obeys some boundedness conditions [52j . 

In short, the result obtained in [53] is that the Cauchy problem for metric f{R) gravity is 
well-formulated and is well-posed in vacuo and with "reasonable" forms of matter (i.e., perfect 
fluids, scalar fields, or the Maxwell field) while for Palatini f{R) gravity, instead, the Cauchy 
problem is not well-formulated nor well-posed due to the presence of higher derivatives of the 
matter fields in the field equations and to the fact that it is impossible to eliminate them |53j . 

Let us consider the scalar-tensor action 



S= d^x 



+ ; (4.23) 



Salgado [50] showed that the corresponding Cauchy problem is well-posed in vacuo and well- 
formulated otherwise. With the exception of w = —3/2, Salgado's results can be extended to 
the more general action 



S= / d^x 



—\l (pVcCp - V{(t)) 



j^gimatter) ^ (4.24) 



which contains the additional coupling function a;((/>) |53j . 

Setting K = 1 in this section and performing a 3 + 1 Arnowitt-Deser-Misner splitting, one 
introduces lapse N , shift A'^*, spatial metric hij, extrinsic curvature Kij, and spatial gradient Qi 
of (/) [21 mi [50]. Assume that a time function t exists such that the spacetime {M,gab) admits 
a foliation with hypersurfaces Sj of constant t with unit timelike normal n". The 3-metric and 
projection operator on St are hab = Qab + naUb and respectively, with 

n^Hc = -1, hab'n!' = habrf = 0, hahbc = hac , (4-25) 

(4.26) 

ds^ = - (A^2 _ Js^^J^.■^ ^^2 _ 2Nidtdx' + hijdx'dx^ (4.27) 
{i,j = 1, 2, 3), while Qc = Dc(j), and the momentum of is 11 = = n'^'Vc(f>- Moreover, 

K,, = - V,n, = + D,N, + D,N,^ , (4.28) 

n = 1 {dt<P + N'Qc) , (4.29) 
dtQi + N'diQi + QAN^ = D, {NH) . (4.30) 



Omitting the calculations, the reduced 3 + 1 equations are [53] 



N 
20 



-(3)ff'j.iV - NKK'j + —6] [2V{(t>) + □(/)] 



dtK + A^'^i^s: +(3) AiV - NKijK'^ 

= — \-2V((P) - 3D(j) + + ^(""^l , (4.32) 
2(p I J 

+ Ij = ^ - 2^(0) + </>y'(</>) + ^ (n^ - Q') ■ (4.33) 

Note that second derivatives of the scalar (j) appear only in the form □(/). For loq = Brans- 
Dicke theory, equivalent to metric f{R) gravity, one can use eq. (I4.33P to eliminate completely 
the d'Alembertian IDcj) from the remaining equations. As a result, the Cauchy problem is well- 
formulated in general and well-posed in vacuo. Further work by Salgado and collaborators |51j 
established the well-posedness of the Cauchy problem for scalar-tensor gravity with a; = 1 in 
the presence of matter, which implies well-posedness for metric f{R) gravity with matter along 
the lines established above. 

Palatini f{R) gravity, instead, is equivalent to an loq = —3/2 Brans-Dicke theory: for this 
value of the Brans-Dicke parameter, the d'Alembertian □(/> disappears from eq. (j4.33p and the 
field 4> is not dynamical (there is no wave equation to govern it): it can be specified arbitrarily on 
a spacetime region provided that its gradient obeys the constraint (I4.33p . Hence, in Palatini f{R) 
gravity it is impossible to eliminate □(/> from the system of differential equations unless, of course, 
□(/> = (including the case of general relativity if (f) = constant). Apart from the impossibility of 
a first-order formulation, one sees that for iv = —3/2 the dynamical wave equation for (p is lost 
completely and this field is non-dynamical. Palatini f{R) gravity has an ill-formulated initial 
value problem even in vacuo and is regarded as a physically unviable theory. 

An alternative approach to the initial value problem is by mapping the equivalent Brans- 
Dicke theory into its Einstein frame representation, in which the (redefined) scalar degree of 
freedom couples minimally to gravity but nonminimally to matter [5l]. This nonminimal cou- 
pling is, of course, absent in vacuo and, from the point of view of the Cauchy problem, plays a 
very minor role in the presence of matter. In this approach, the non-dynamical role of the scalar 
is even more obvious, and the conclusions above are reached by using well-known theorems on 
the initial value problem of general relativity with a scalar field |54j . 

The problem with Palatini f{R) gravity has been noticed with an entirely different approach, 
i.e., matching static interior and exterior solutions with spherical symmetry [56j (other problems 
are reported in Refs. [71 \57\ [58]). 

The field equations of Palatini f{R) gravity are second order PDEs in the metric. Because 
/ is a function of R, which is an algebraic function of T due to eq. (j2.12p . the right hand side 
of eq. ()2.15p includes second derivatives of T. But T contains derivatives of the matter fields 
up to first order, hence eq. (|2.15p contains derivatives of the matter fields up to third order. 



This is in contrast with the situation of general relativity and most of its extensions, in which 
the field equations contain only first order derivatives of the matter fields. As a consequence, 
in these theories the metric is generated by an integral over the matter sources and, therefore, 
discontinuities (or singularities) in the matter fields and their derivatives do not imply unphysical 
discontinuities of the metric. In Palatini f{R) gravity, instead, the algebraic dependence of the 
metric on the matter fields creates unacceptable discontinuities in the metric and singularities in 
the curvature, which is what is found in [56]. So, both the failure of the initial value problem and 
the occurrence of curvature singularities in the presence of discontinuities in the matter fields 
or their derivatives can be traced to the fact that the scalar degree of freedom is non-dynamical 
and is related algebraically to T. A possible cure is to modify the gravitational sector of the 
action to raise the order of the field equations. 

5 Conclusions 

We are now ready to summarize the situation of f{R) gravity. Let us stress once again that we 
regard these theories more as toy models, and as proofs of principle that modified gravity can 
explain the observed acceleration of the universe without dark energy, than definitive theories. 

• Metric f{R) gravity: models exist that pass all the observational and theoretical con- 
straints. An example is the Starobinsky model [59] 



All the viable models require the chameleon mechanism in order to pass the weak-field 
limit tests. A condition that must be satisfied by all metric f{R) theories in order to 
avoid the Dolgov-Kawasaki local instability is f"{R) > 0. The condition ()4.12p must be 
satisfied for the stability of a de Sitter space. The biggest problem is whether curvature 
singularities exist for relativistic strong field stars. 

• Palatini f{R) gravity: these theories suffered multiple deaths; they contain a non- 
dynamical scalar field, the Cauchy problem is ill-posed, and discontinuities in the matter 
distribution generate curvature singularities. 

• Metric-afRne gravity: this class of theories is not yet sufficiently developed to assess 
whether it is viable according to the criteria listed here, and its cosmological consequences 
are unexplored. 

In conclusion, f{R) theories have helped our understanding of the peculiarities of general 
relativity in the broader spectrum of relativistic theories of gravity, and have taught us about 
important aspects of its simple generalizations. They even constitute viable alternatives to dark 
energy models in explaining the cosmic acceleration, although at present there is no definite 
prediction that sets them apart once and for all from dark energy and other models. 
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